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EXPERIMENTAL AND THEORETICAL DETERMINATION 
OF THEPMAL STRESSES IN A FLAT PLATE 
By Richard R. Heldenfels and William M. Roberts 

SUMMARY 


Theimal stresses induced in a flat, rectangular, 75S-T6 aliaminum- 
alloy plate by nonuniform heating are determined both experimentally and 
theoretically. The characteristics of commercially available bonded 
resistance wire strain gages are first investigated to determine their 
suitability for measuring stresses under simple conditions of stress 
and temperature. The gages are then used to measure thermal stresses 
in the flat plate in order to study their suitability under more com- 
plicated conditions. The experimental results are found to be in 
satisfactory agreement (within percent of the maximum calculated 
stress) with an approximate theoretical solution of the problem. 


INTRODUCTION 


Aerodynamic heating of aircraft flying at supersonic speeds has 
become important in structural design because material properties are 
changed by exposure to elevated temperatures, and thermal stresses which 
result from nonuniform heating may have a significant effect on the 
strength of the structure. Methods for determining thermal stresses 
are therefore required if an efficient design is to be attained. 

References 1 and 2 considered the thermal-stress problem in air- 
craft structures theoretically and presented some procedures for thermal- 
stress analysis. The methods presented varied in their accuracy and 
involved various assumptions and approximations. These theoretical 
approaches alone, therefore, are not sufficient, and experimental methods 
are needed to check the accuracy of approximate thermal-stress -analysis 
procedures and to help analyze those structures which are too complex 
for theoretical analysis. 

An investigation has therefore been made of the characteristics 
of commerically available bonded resistance wire strain gages to deter- 
mine their suitability for measuring stresses under simple conditions of 
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stress and temperature, and the results are presented in this paper. 

The strain gages are used to measure the thermal stresses induced in a 
flat plate subject to a known, nonuniform temperature distribution to 
study their suitability under more complicated conditions. The stresses 
so measured are then compared with an approximate theoretical solution, 
the details of which are included in the appendix. 


SYMBOLS 


a 

b 

E 

Eo 

T 

To 


U 

X, y 

a 


€ 

€ * 
X 



half-length of plate, in. 
half -width of plate, in. 
modulus of elasticity, psi 

modulus of elasticity at room temperature, psi 
temperature, 

temperature at which thermal stresses are assumed to be zero. 

Op 

temperature difference between center line and longitudinal 
edge of plate, ^F 

strain energy, in. -lb 

coordinate axes 

coefficient of thermal expansion, — 

in. -^F 

total- strain reading of gage, in. /in. 
net strain, € - X, in. /in. 

zero shift of strain gage due to temperature changes, in. /in. 
Poisson ^s ratio 

longitudinal direct stress, psi 
transverse direct stress, psi 
shear stress, psi 
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<I) stress function 




differential operator, 


differential operator. 


^ + 2 


STRAIN-GAGE CHARACTERISTICS 


An ideal strain gage, when mounted on a test specimen, would produce 
stable, accurate, and reproducible readings under all conditions of 
stresses and temperatures. That is, the reading would not drift under 
constant stress and temperature, the reading would be unique for each 
combination of stress and temperature, and the reading would be the same 
for different gages of the same type. If such a gage were available, 
the measured strain might be simply related to a uniaxial stress in the 
elastic range as follows: 


a = = E(€ - X) 


where 

a stress, psi 

E modulus of elasticity of the specimen which is a function of 

temperature, psi 

e* net strain, in. /in. 

€ measured strain, in. /in. 

X zero shift due to temperature change, in. /in. (called zero 

shift because it is the measured strain due to temperature 
changes under zero stress conditions) 

If bonded resistance wire strain gages could closely approximate such 
behavior, they would provide a simple means of measuring thermal stresses. 
A study was therefore made to determine the characteristics of one type 
of wire strain gage and its suitability for measuring thermal stresses. 

Commercially available Baldwin SR-^ type AB-7 strain gages were 
selected for investigation. These -1 - inch-gage-length bakelite gages 
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were chosen because of their good elevated-temperature characteristics 
(ref. 3 ) and because subsequent tests required a short gage length. The 
AB-7 gage was selected in preference to self -compensating gages because 
the latter were not available in sufficient quantity when this study 
was undertaken. The higher cost of the self -compensating gages also 
discourages their use when large numbers of gages such as were used 
herein are needed. 

The characteristics of the AB-7 gages were studied at temperatures 
between 80° F and 300^ F and at various stress levels in the elastic 
range of 75S-To aluminum alloy. The study was carried out by mounting 
a gage on each side of a 75S-T6 aluminum-alloy block. The gages were 
mounted in accordance with procedures recommended by the manufacturer 
and were read separately on a type K Baldwin SR-4 strain indicator . A 
pure compression load was applied to the block by means of a hydraulic 
testing machine and temperatures were controlled by an electric furnace 
(ref. 4). Throughout these and all subsequent tests, a consistent 
wiring arrangement was used so that wiring effects would be included in 
the calibration. 

The specimen and the mounted gages were first exposed to conditions 
of constant stress and temperature to investigate the stability of the 
gages — that is, to find whether the gage reading changed with time and, 
if so, how much it changed. Changes of 10 microinches per inch, probably 
due to creep of the gage and specimen, were detected for exposures of 
approximately l/2 hour at 300^ F and 12 ksi but were considered small 
enough to be neglected for this condition and for longer times at lower 
stresses and temperatures. The gages are thus deemed sufficiently stable 
for short-time tests. 

Accuracy and reproducibility were next investigated by obtaining 
compressive stress-strain curves in the elastic range for teii 5 )eratures 
up to 300 ^ F. Figure 1 shows typical stress -strain curves at three 
temperatures for an AB-7 mounted on 75S-T6 aluminum alloy. The 

curves of figure 1 were obtained by measurement in the compression range 
and extrapolation into the tension range by assuming that the modulus of 
elasticity is the same for tension and compression. The curves of fig- 
ure 1 are defined by two quantities: the slope of each line, which is 
actually the modulus of elasticity of the specimen, and the spacing along 
the horizontal axis, which is the zero shift. 

The ratio of the nodulus of elasticity of 75S-T6 aluminum alloy at 
the test temperature to its modulus at room temperature, both measured 
by the AB-7 gages, is plotted in figure 2. Included in the figure are 
the average values and scatter band of the modulus ratio of 75 S-T 6 alumi- 
num alloy as determined by a different method and reported in reference 4. 
The modulus ratios obtained by means of AB-7 gages are seen to be within 
the accuracy of those of reference 4 and the gage factor is apparently 
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constant. Reference 3 also found that the gage factor of AB-7 gages was 
constant in this temperature range. Reproducibility was indicated since 
different gages gave substantially the same modulus ratios. 

The zero shift differed from gage to gage on the same specimen and 
was dependent upon the previous temperature history of the gage. The 
response of a particular gage to temperature changes was found, however , 
to become fairly consistent if the gages were mounted, cured at a tem- 
perature above that at which they were to be used, and then put through 
several temperature cycles. Even after this treatment, however, the 
responses varied from gage to gage within the range shown in figure 3- 
The data of figure 3 were obtained from 70 of the 72 different gages 
used in the investigation, the zero shift of the other two being dis- 
carded because of erratic behavior. The largest deviation of the 
AB-7 gages from ideal performance then was in the gage-to-gage variation 
of the zero shift which was apparently dependent on small variations in 
gage construction and installation. 

These tests were concerned with uniaxial stresses . If more com- 
plicated states of stress are of interest, the measuring gages must be 
arranged in rosettes. In such cases, the stresses can be calculated by 
the applicable rosette formula if the net strain is used for each arm 
of the rosette. Reference 5 presents detailed information on the use of 
rosette gages and shows that the transverse sensitivity is a small 
quantity. A confirmatory investigation at room temperature indicated 
that the transverse sensitivity of AB-7 strain gages could be neglected. 

The AB-7 wire strain gages therefore have satisfactory character- 
istics for the measurement of stresses under conditions of uniform stress 
and temperature, provided that the strain gages are properly mounted on 
the specimen and are cured and subjected to several temperature cycles 
and provided that the variation of the zero shift with temperature is 
determined for each individual gage. A complete stress-strain curve of 
the type shown in figure 1 is not necessary for each gage since the 
variation of elastic modulus with temperature can be obtained from any 
suitable source. In many tests, calibration of all of the gages mounted 
on the specimen may be impractical, and in such cases, an average zero 
shift can be used if accuracy can be sacrificed. The present investiga- 
tion shows, however, that the use of an average zero shift would have 
given errors up to 1 ksi. 


MEASUREMENT OF THERMAL STRESSES IN THE TEST PANEL 

The tests described in the preceding section showed that the 
AB-7 gages had satisfactory characteristics for the measurement of 
thermal stresses under simple states of temperature and stress. The 


6 


NACA TN 2769 


suitability of the gages for the measurement of thermal stresses under 
more complicated conditions was therefore investigated by means of a 
test panel designed to provide a two-dimensional distribution of thermal 
stress under steady-state conditions. A simple configuration, arranged 
for experimental convenience, was chosen so that the accuracy of the 
experimental measurements could be verified by comparison with a theo- 
retical solution. This configuration was also suitable for an experi- 
mental and theoretical investigation of thermal buckling of plates 
(ref. 6). 

A flat, rectangular plate of ^ - inch- thick T5S-T6 alviminum alloy, 

measuring 24 by 36 inches, to which apparatus was attached for inducing 
a steady-state temperature gradient and measuring the resultant tempera- 
tures and stresses (figs. 4 and 5) vas used. The temperature gradient 
was induced in the plate by heating its longitudinal center line and 
cooling its longitudinal edges. The heat was supplied by Nichrome wire 
imbedded in porcelain cement and was removed by water flowing through 
two plastic tubes glued to the plate. The size and material of the 
attached apparatus were such that negligible resistance was offered to 
thermal expansion of the plate. The sides and ends of the plate were 
insulated with asbestos (see fig. 4) to minimize heat loss through these 
surfaces and to insure that the instrumentation was at the same tempera- 
ture as the plate. The flow of cooling water was controlled so that the 
temperature distribution was symmetrical about the longitudinal center 
line and essentially constant in the longitudinal direction. No measur- 
able temperature variation occurred through the thickness of the plate. 

Since the temperature distribution, and hence the stress distribution, 
was symmetrical about both center lines of the plate, strain gages were 
mounted on only one quadrant of the panel. At each point where strain 
measurements were taken, three AB-7 strain gages were arranged in a fan 
rosette and a similar group of strain gages was mounted on the reverse 
side. Iron-coiistantan thermocouples were mounted at the center of each 
fan on one side of the plate. 

The thermal stresses induced in the plate tended to deform it out 
of its plane, a combined bending -buckling phenomenon associated with the 
initial imperfections of the plate. Since the present tests were con- 
cerned with plane stresses, this deformation was minimized by mounting 
the plate in a jig designed to provide the required restraint without 
introducing extraneous stresses. To further insure that the experimental 
stresses would closely approximate a state of plane stress, the readings 
of strain gages on opposite sides of the plate were averaged before the 
stresses were computed. 

Three tests were made and in each the panel was heated until the 
temperature distribution stabilized; then the strain gages were read. 

The strain and temperature data were processed as follows: 
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( 1) The temperature distribution was plotted and the average 
temperature of each strain gage determined. 

(2) The zero shift of each gage was determined from its calibration 
curve. 

(3) The net strain of each gage was obtained by subtracting its 
zero shift from its total strain. 

(4) The net strains were plotted and the corrected net strains at 
the common point of each rosette were interpolated. This correction 
was necessary to obtain accurate results because a large strain gradient 
existed over some areas of the panel and the three gages of the rosette 
were necessarily located some distance apart. 

(5) The corrected net strains from gages on opposite sides of the 
plate were averaged to get the midplane strain. 

(6) The three stresses at each point (direct stresses in longi- 
tudinal and transverse directions and the shear stress) were computed 
from the three corrected net strains and the known properties of the 
material. 

(7) The resulting stresses were then adjusted to a temperature 
difference of 150^ F between the center line and longitudinal edges of 
the plate by assuming that the stresses were a linear function of the 
temperature difference. This adjustment was made to provide a common 
basis for comparing the three tests which were made at different tempera- 
ture differences. Actually, the stresses do not vary linearly with 
temperature difference because of the changes in material properties; 
however, in the temperature range of the tests , the departure from 
linearity was very small because the decrease in the modulus of elasticity 
with temperature increase was offset by the increase in the coefficient 

of thermal expansion with temperature increase. 

The results are shown in figures 6 to 9^ in which the experimental 
results are compared with the theoretical results described in the next 
section. Figure 6 shows the transverse temperature distribution obtained 
at all longitudinal stations. Figure 7(a) shows the transverse distribu- 
tion of the longitudinal direct stress at each of the transverse stations, 
and figure 7(b) shows the longitudinal distribution of the longitudinal 
direct stress at each of the longitudinal stations. Similarly, the shear 
stresses are shown in figures 8(a) and 8(b) and the transverse direct 
stresses are plotted in figures 9(a) and 9(b). The experimental points 
shown are the average of the three tests; separate points are not shown 
for each test since in most cases they are too close together to be 
easily distinguishable. 
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CALCULATION OF THERMAL STRESSES IN THE TEST PANEL 


The stress distribution in the test panel was calculated by applica- 
tion of the theory of elasticity to provide a theoretical check on the 
measured stresses. The analysis is simplified by assuming that the test 
panel is a thin, flat, rectangular plate which is in a state of plane 
stress and which does not buckle out of its original plane and that the 
properties of the material are constant in the temperature range of 
interest . 

The plate is free to expand in all directions; no external loads 
axe applied; and all internal stresses are caused by a nonuniform tem- 
perature distribution. Pertinent dimensions and the coordinate system 
of such a plate are shown in figure 10 and the temperature distribution 
is shown in figures 6 and 11. 

An exact solution of this problem can be obtained in the form of 
an infinite series. Such a solution, however, is inconvenient to 
evaluate numerically, and an approximate solution which is developed in 
the appendix is sufficiently accurate for general verification of the 
test results. Such an approximate solution is based on the assumption 
that the variation of the longitudinal direct stress Ox in the trans- 
verse direction is the same as that given by the analysis of a very long 
plate in reference 7. The transverse distribution of shear stress r^y 
and the transverse distribution of direct stress Oy are determined 
from the longitudinal direct stresses a^, equilibrium equations, and the 

stress-free boundary conditions on the longitudinal edges. The longi- 
tudinal variations of these stresses are obtained by application of the 
principle of minimum complementary energy. 

The dimensions, temperature difference, and material properties used 

are 


a = 18 in. 
b = 12 in. 
Ti = 150 ° F 


E = 10.4 X 10 


6 


psi 


a = 12.7 X 10 


-6 


in. 

in.-°F 
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and the results of the theoretical analysis of the test p?nel are 
CT^ = 1651(7 - 6) E - 0.0987 sinh(0.1759x)sin(0.0923x) - 


0.1528 cosh(0.1759x)cos(0.0923x 


>] 


Ty.y = 82.55(7^ - 12y) 1^.0326 cosh(0.1759x)sin(0.0923x) + 
0.3600 sinh(0.1759x)cos(0.0923xj 

Oy = 2.752(7^ - l8y2 + 864 )[o .2750 sinh(0.1759x)sin(0.0923x) - 


0 .6634 cosh( 0 . 1759x ) cos ( 0 . O923X ) 


(0 X 18; 

0 ^ y ^ 12) 


These equations give the stresses in the first quadrant; the stresses 
in the other quadrants can he obtained by noting that the direct stresses 
cr^ and Cy are symmetrical and the shear stress is antisymmetrical 

about both the longitudinal and transverse center lines • 


RESULTS AND DISCUSSION 
Temperature Distribution 


The temperature increase of the instrumented quadrant of the test 
panel is shown in figure 6. The transverse distribution is shown for 
only one station; the others are not shown since they are identical. 

The temperature of the panel when the test began (80*^ F) was the tempera- 
ture at which the initial gage readings were taken and thus is the tem- 
perature at which the thermal stresses were assumed to be zero. 

The straight line drawn through the test points is the temperature 
distribution used to calculate the theoretical theimal stresses and 
corresponds to simple conduction between a concentrated heat source along 
the longitudinal center line and a concentrated heat sink along each 
longitudinal edge. The temperature difference shown between the center 
line and the edge is 150° F. The actual temperature distribution, given 
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by the test points, varied slightly from a straight line because of the 
finite area of the source and sink and because of heat loss through the 
sides and ends of the plate, but the measured points are within a few 
degrees of the straight line. 


Longitudinal Direct Stresses 

The longitudinal direct stresses induced in the plate by the tem- 
perature distribution of figure 6 are shown in figure 7- figure 7(9-) 

the transverse distribution of calculated and measured stresses is shown 
for four stations, and in figure 7(b) the same data are replotted to show 
the longitudinal distribution at three stations. The stresses are shown 
for only one quadrant of the plate, but the complete distribution is 
symmetrical about both center lines. 

The results show that the hot parts of the plate near the longi- 
tudinal center line (y = O) are in compression and the cold parts near 
the longitudinal edges (y = 12) are in tension, as would be expected, 
since the cold parts restrain the greater expansion of the hot parts. 

The theoretical results shown are seen to be in good agreement with the 
test data. 


Shear Stresses . 

The shear stresses induced in the plate are shown in figure 8, in 
which the presentation is similar to that in figure 7* Again the stresses 
are shown for only one quadrant of the plate, but in this case the com- 
plete distribution is antisymmetrical about both center lines if a 
uniform sign convention is used for shear stresses. The agreement 
between calculated and measured shear stresses is fairly good but not 
so good as for the longitudinal direct stress. 


Transverse Direct Stresses Cy 

The transverse direct stresses induced in the plate are shown in 
figure 9 w’here the presentation is similar to that in figures 7 8* 

Results are shown for only one quadrant, but the complete distribution 
is symmetrical about both center lines . The agreement between calculated 
and measured values again is not so good as for the longitudinal direct 
stress . 
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Agreement Between Theory and Experiment 

The general agreement between the measured and computed stresses is 
satisfactory with the difference being less than l /2 ksi for all test 
points except one. Both the test data and theory are subject to inac- 
curacies, the magnitudes of which are not precisely known. The test 
and theory, however, mutually support one another and are of about the 
same order of accuracy, with the measurements probably being more 
accurate than the theoretical values. The fact that good agreement exists 
between calculated and measured longitudinal direct stresses and only 
fair agreement between calculated and measured transverse direct and 
shear stresses can be attributed to the errors introduced by the approxi- 
mate analysis rather than errors in measurements. 

A long thin plate subjected to a temperature distribution that is 
symmetrical about the longitudinal, center line and independent of the 
distance along the plate length is analyzed in reference 7 and is shown 
to have induced longitudinal direct stresses but no transverse 

direct stresses or shear stresses ^xy ^ short plate, such as 

investigated in the present paper, the condition that the edge stresses 
vanish necessitates the modification of the longitudinal direct stresses 
and introduces transverse direct stresses and shear stresses into the 
problem. 

In the approximate -analysis of the present paper the transverse 
variation of the longitudinal direct stress given by the long-plate 
solution is assumed to be valid for short plates and a correction is 
determined for the longitudinal variation of these stresses by the 
principle of minimum complementary energy. Since the long-plate solu- 
tion does not involve transverse direct stresses and shear stresses, the 
corrected longitudinal direct stresses can be expected to be nearer the 
true values than the transverse direct stress Oy and the shear stress 

*^xy vhich are themselves corrections. Additional accuracy could be 
obtained by making additional corrections; however, the analysis used is 
estimated to be within 1/2 ksi ( ^5 percent of the maximum calculated 
stress) of the exact solution in the region where the stresses were 
measured. Since this is of the same order of magnitude as the experi- 
mental accuracy, a more refined theoretical approach is not required. 


CONCLUDING REMAEKS 


A technique has been described for measuring thermal stresses with 
commercially available bonded resistance wire strain gages. In such 
applications, wire strain gages cannot be expected to yield as accurate 
results as in more conventional uses because the temperature changes are 


12 


NACA TN 2769 


large and there is no simple way to compensate for the many errors thus 
introduced. The largest deviation from ideal performance was in the 
gage-to-gage variation of the zero shift due to temperature changes which 
was apparently dependent on small variations in gage construction and 
installation. This difficulty was overcome in the technique described 
by measuring the zero shift of each gage after it was installed on the 
test specimen and cured by several cycles of temperature higher than 
that reached during the test. Such a procedure may be impractical in 
many cases, since it involves uniform heating of the complete test 
specimen over the entire temperature range of interest. If this calibra- 
tion is impractical and accuracy will permit, an average value of the 
zero shift obtained from a separate test may be used. In this case, 
however, errors up to 1 ksi may be expected in the range of conditions 
investigated in this paper. 

The technique described was applied to the measurements of thermal 
stresses induced in a flat plate subjected to a known steady-state 
temperature gradient. The stresses so measured were compared with an 
approximate theoretical solution and were found to be in satisfactory 
agreement, that is, within about ±5 percent of the maximum calculated 
stress . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , June 9 , 1952 
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APPENDIX 

DERIVATION OF APPROXIMATE FORMULAS FOR 
THERMAL- STRESS DISTRIBUTION 


The rectangular plate shown in figure 10 is analyzed under the 
following conditions: 

(1) The temperature distribution is 'unifonn over the plate thickness 
and varies over the plate surface in the following manner: 


T = Tq + XY 


(Al) 


where X is a function of x only and Y is a function of y only. 
The uniform temperature distribution Tq is superimposed on the plate 
without affecting the problem — that is, a uniform temperature distribu- 
tion produces no thermal stresses because of the next condition. 

(2) No external loads are applied and the plate is free to expand 
in all directions so that all stresses are due to the nonuniform temper- 
ature distribution. 


The following assumptions are also made: 

(1) The plate is thin and may be considered to be in a state of 
plane stress. 

(2) The properties of the material do not vary with temperature 
and the stresses are within the elastic range. 


It is convenient to introduce the Airy stress function 0 such 

that 


cr 


X 



" ■ ^x Sy 





(A 2 ) 
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The exact solution of the stated problem is then given by the function ^ 
that satisfies the following differential equation (ref. 7)- 


= -EaV^T 


(A3) 


A solution of equation (a 3 ) can be obtained in the form of an infinite 
series; however, such a solution is inconvenient to evaluate numerically. 
The type of series encountered may be visualized by noting that 0 is 
mathematically equivalent to the deflection of a rectangular plate of 
uniform thickness with clamped edges and under a lateral load propor- 
tional to V^T. 

It is further assumed that, as a first approximation, the Airy 
stress function 0 can be expressed as 


0 = f g 


(aU) 


where f is a function of x only and 
Then 


is a function of y only. 

•N 


= fg" 

T^y = -fg’ 
CTy = f’g 


(A5) 




The primes indicate differentiation of the fvinction (f or g) with 
respect to its assigned variable (x or y) . 


The problem is solved by selecting an appropriate function g and 
then using the principle of minimum complimentary energy to determine a 
function f that gives the best approximation to the exact solution of 
equation (A3). The degree of approximation involved is dependent upon 
the function g selected. A good selection is a function proportional 
to the stress function of an infinitely long plate (a = “) subject to 
the temperature distribution T = Y. In this case, the stress distribu- 
tion is independent of x; therefore, equation (A3) requires that 


,IV 


Y" 


(A6a) 
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or 


g" oc Y + Bjy + Bg (A 6 b) 


Since, from equation (A5), is proportional to g", the distribution 

of is proportional to the temperature distribution Y plus a linear 

distribution + B 2 required to maintain thrust and moment equilib- 

rium on the cross section. These results are the same as those obtained 
for an infinitely long plate by a somewhat different process in refer- 
ence 7- The function g is obtained by integrating g” twice and the 
two additional constants of integration are determined from conditions 
at the edges y = ±b. This process is illustrated subsequently by a 
particular example. In the following paragraphs, the equations for 
determining f are developed with the assumption that g is some known 
function of y. 

The complementary energy of a heated plate may be written as 






- 2 f-io’^o’y + 2(1 + 


2EoT(a^ 



dx dy 


(AT) 


where the term involving oT expresses the strain energy resulting 
from thermal expansion. If equations (Al) and (A5) are substituted 
into equation (AT), the complement ary- energy expression becomes 



A2f 


m 2 


- anA^ff" + 2(1 + ii)A4f'^ + 




(A 8 ) 
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The complementary energy is now minimized hy setting the first 
variation of U, obtained by the usual procedures of the calculus of 
variables, equal to zero. Since g is a known function of y, the 
result, after integration by parts, is 
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5U 


' I J* ri*' - 2 [^1. + ^‘(A 3 + 

Ea(Aj + AgX + A3X”)| Sf dx + 

I kiA^ + 2(1 + n)A]^f - Agf "' - EoAsX' 

||jf2^" ■ Ea(AY + AsXyjsfJj 


(A9) 


The terms under the integral sign yield the differential equation 
for the function f 


A2f 


IV 




fM 


+ A;^f = -Ea^A:^ 


AgX + A3X") (AlO) 


The remaining terms set equal to zero are the boundary conditions. The 
problem as stated requires that the edges of the plate be stress -free, 
so that physically the boundary conditions are: 

At X = la, 


c^x = ■’’xy " ^ (Alla) 

and at y = lb, 

Cy = = 0 (Allh) 

Equations (Alla) require that f and f* vanish at x = la, so that 
the requirements of the boundary terms of equation (A^) are satisfied. 
Equations (Allh) give boundary conditions which must be considered in 
the selection of g and require that, from equation (a8), 


(A 12 a) 


^3 ” 

A^ = 0 


(A 12 b) 
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Equation (A12a) shows that the approximate solution is not dependent 
on Poisson's ratio p, a result that is not surprising since the exact 
solution does not involve p. The differential equation for f then 
"becomes 


Agf^^ - 2\f" + Ajf = -Ea^AgX + A 3 X") 


(AI3) 


Equation (AI3) is now solved for the temperature distribution shown 
in figure 11. This distribution, which does not vary with x, has a 
discontinuity along the x-eixis and cannot be defined by a single, simple 
analytical expression. The stress distribution in the plate due to this 
temperature distribution is symmetrical about the plate center lines. 

The analysis can therefore be confined to the first quadrant of the 
plate. The temperature distribution in this quadrant is then given by 


X = T^ = Constant 

Y= 1 -g 

T = To + Ti(l - 2) 


(Al4) 


The fionction g is found from the criterion of equation (A6b) and 
the requirement that the stress distribution is symmetrical about the 
plate center line; that is, 


g" oc Y + B2 = 1 - + Bg 


(AI5) 


or, since an arbitrary constant can be removed without affecting the 
solution, 


g 





(Al6a) 


Then 


g’ 





(Al6b) 
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and 


g 




( A16c) 


The boundary conditions of equation (Allb) should determine the 
constants of integration in equations (A16); however, only one quadrant 
is being analyzed and the conditions at y = -b are not applicable. 
Conditions of symmetry are used instead — that is, and Oy are 

made symmetrical about the x-6Lxis. The boundary conditions required to 
define g are then 


g'(0) = 0 

g(b) = 0 
g'(b) = 0 




J 


(AIT) 


in which case 

(A18) 

Since Y and g are given by equations (Al^l) and (A18), the 
integrals k-^ to Aq can be evaluated to determine the coefficients 

of equation (AI3). After some simplification the result is 



13b^f^^ - 84b^f" + 420f = 420EaTi 


(AI9) 


The solution of equation (A19) can be written as follows: 


f = EaT^^l + C]^ sinh k^x sin k2X + C2 sinh k2^x cos k2X + 


cosh k]^x sin k2X + cosh k^x cos k2X^ 


(A20) 
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where 


k 


2 

1 


21 4- Vi36^ 
13b2 


and 



-21 -I- ^1365 
13b^ 


The boundary conditions of equation (Alla) determine the constants 
of integration, but conditions of symmetry may be used instead of the 
conditions at x = -a since the analysis is confined to the first 
quadrant. The required relations are then: 


f'(0) = 0 
f (0) = 0 
f(a) =0 
f'(a) = 0 




> 


(A21) 


The substitution of equation (A20) into equations (A2l) gives four 
equations for the constants of integration to which may be 

solved to yield -v 

sinh k^a cos - k^ cosh k^a sin k^a 

Cn 

k^^ sin k 2 a cos k 2 & + k 2 sinh k^a cosh k^^a 


C2 



(A22) 




cosh kj^a sin k 2 & + kg sinh k^a cos kgS 
k]^ sin kga cos kga + kg sinh k-^a cosh k^a 


The following expressions for the stresses can be obtained by 
substituting equations (A18) and (A20) into equations (A5): 
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EoT-i 

= 

X 2 


^2 ^ sinh sin k2- 


X + 


Cj^ cosh kj^x cos k2^^ 


EaT] 
’’xy ~ 


>[(lf - i(' 


sinh cos + 


(A23) 


Cg cosh k^x sin k2^^ 


EaT^ ^2 

y " 2 6 


2(1 


, w '- 3 gT -‘ 


sinh k^x sin k2X + 


C0 cosh k^x cos k2^^ 


where 




C6 = C]_k3^ - C4k2 

C7 = C^(k^^ - k2^) - 2C4kik2 

Cq = 2Cikik2 + Ci^^k^2 _ 
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Figure 1.- Typical stress-strain curves for wire strain gages mounted on 

75S-T6 alumimmi alloy. 
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Figure 2 
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- Ratio of modulixs of elasticity at elevated temperature to that 
at room temperature for 75S-T6 alximinum alloy. 
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Figure 3 .- Zero shift, due to temperature change, of TO wire strain gages 

moimted on 75 S-T 6 aluminum alloy. 
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Figure k.- Sketch of thermal-stress test specimen. 
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Figiire 5 *- Photograph of test setup. 
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Figure 6 . - Transverse temperature distribution. 
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(a) Transverse distribution. 

Figure 7.- Longitudinal direct stresses induced in test specimen by the 
temperature distribution of figure 6. 
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(b) Longitudinal distribution. 
Figure 7* - Concluded. 
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Figure 8 
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(a) Transverse distribution. 

. - Shear stresses induced in test specimen by the temperature 
distribution of figure 6 . 
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(b) Longitudinal distribution. 
Figure 8 .- Concluded. 
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(a) Transverse distribution. 

Figure 9.- Transverse direct stresses induced in test specimen by the 
temperature distribution of figure 6 . 
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(b) Longitudinal distribution. 
Figure 9. - Concluded. 
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Figure 10.- Dimensions and coordinate system of plate analyzed in the 

appendix. 




Figure 11.- Transverse temperature distribution on plate analyzed in the 

appendix. 
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